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Abstract 

General structures of the charged lepton and the neutrino mixing matrices leading to tri- 
bimaximal leptonic mixing are determined. These are then integrated into an 50(10) model 
within which detailed fits to fermion masses and mixing angles are given. It is shown that one 
can obtain excellent fits to all the fermion masses and quark mixing angles keeping tri-bimaximal 
leptonic mixing intact. Different perturbations to the basic structure are considered and those 
which can or which cannot account for the recent T2K and MINOS results on the reactor mixing 
angle are identified. 
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I. INTRODUCTION 



The Tri-bimaximal (TBM) leptonic mixing [T] provides an important clue in search of 
possible flavour structure [2] governing the leptonic masses and mixing angles. It predicts 
sin^^5^2 = 1/3 and sin^^23 — 1/2 respectively for the solar and the atmospheric mixing 
angles both of which agree nearly within la with the latest global analysis [31 H] of the 
neutrino oscillation data. TBM also predicts vanishing reactor mixing angle ^^3. This can 
be reconciled with the latest T2K [5J (MINOS f6]) results at 2.5a (I.60-) and with the 
global analysis [3l H] at about 3a. This suggests that the TBM may be a good zeroth 
order approximation which needs perturbations affecting mainly the reactor mixing angle 
^^3. While such perturbations may arise from some underlying flavor symmetry, see [7j for 
examples, it would be more appropriate if some independent mechanism like grand unifled 
theory (GUT) [S] governs these perturbations. We wish to analyze here the TBM structure 
and perturbations to it within a grand unified model based on SO{10) gauge symmetry. 

Incorporating the tri-bimaximal mixing into GUTs, particularly based on the 5*0(10) 
gauge group is quite challenging, see [9| for some examples. Since all fermions in a given 
generation are unified into a single 16 dimensional irreducible representation of 5*0(10), 
imposition of the TBM structure on the leptonic mass matrices also constrains the quark 
mass matrices. It is not clear if the requirement of the exact tri-bimaximal mixing among 
leptons would be consistent with a precise description of the quark masses and mixing. We 
suggest a general method of incorporating the exact TBM structure within 50(10) and use 
it to obtain quantitative description of the fermion masses and mixing. 

Our approach is purely phenomenological. We do not use any fiavour symmetry but 
determine the most general structure of the leptonic mixing matrices required for obtaining 
the tri-bimaximal mixing. We then try to integrate this structure into 5*0(10) and discuss 
numerical fits to the quark masses and mixing leaving the TBM intact. Then we discuss 
perturbations to the basic structure and their effects on observables. 

II. LEPTONIC MIXING MATRICES AND TBM 

We shall derive general forms for the neutrino mass matrix and the left handed 
charged lepton mixing matrix Ui which lead in the fiavour basis to a neutrino mass matrix 
A4uf exhibiting the TBM structure. We define the TBM structure for Ad^f as follows: 



. / 2/1 + /2 /2 - /l /2 - /l \ 

3 /2-/1 l(/l + 2/2 + 3/3) l(/l + 2/2-3/3) 

\f2-f1 Ufi + 2/2 - 3/3) i(/i + 2/2 + 3/3) J 



(1) 



where /i,2,3 are complex neutrino masses. This matrix is diagonalized by 



UpMNS — OtbmQ 5 



(2) 
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where Q is a diagonal phase matrix and 
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A more general definition of TBM structure would be to replace Aiuf and Upmns above 
by PiM.i,fPi and P^Upmns, where Pi denotes a diagonal phase matrix. Since Pi can be 
rotated away by redefining the charged lepton fields, we shall refer to TBM structure as the 
one defined by Eqs. ( [lp| . 

It is known [TD] that Ai uj in Eq. ([T]) is invariant under a Z2 x Z2 symmetry. The elements 
of the Z2 X Z2 are defined as 
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and satisfy 



(4) 



(5) 



Si3M.fS2,3 = M.f 

We shall exploit this Z2 x Z2 invariance in arriving at the structure of Ui. As noted in 
[TT] . one can always choose a specific basis in which exhibits the TBM structure and is 
thus invariant under Z2 x Z2: 

Sl,M,S2,3 = (6) 

If in an arbitrary basis is not invariant under Z2 x Z2 then one can go to a new basis 
with Ml, = U'^MlU and choose U in such a way that 



with a diagonal matrix with real positive elements and 



U,. = O 



TBM 



P . 



(7) 



(8) 



P being a general diagonal phase matrix. Let Ui denote the mixing matrix among the left 
handed charged leptons in a basis in which M^y is Z2 x Z2 symmetric. If such a defined Ui 
itself is Z2 X Z2 symmetric, i.e. satisfies 



SjnUiS2,3 — Ui 



(9) 



then J^uf will also satisfy Eq. ^ and thus would exhibit the TBM structure of Eq. ([T]). 
This follows trivially from the definition 



M 



UfM^Ui 



(10) 



after using Eqs. Thus the Z2 X Z2 invariance of Ui is sufficient to ensure the TBM for 
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Above equation allows us to determine TBM preserving class of Ui in a basis with M^, 
satisfying Eq. ([6]). 5*3 invariance corresponds to imposing the fi-r interchange symmetry on 
Ui. The 5*2 invariance further requires [//^ = Ui and that the sum of elements in each of its 
raw must be equal. Such a Ui can be parameterized as 



Ui 



(11) 



where Pi = diag.(l, e*'^, e*'^) is a diagonal phase matrix and 
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(12) 



with tan 6' = — 2v^cos/3. Ui is thus fully determined by three phase angles a, (3 and 5. 

The form of Ui as given above is the most general one required in order to obtain TBM in 
the basis with M^, satisfying Eq. ([6]). The generality is proved by noticing that the Z2 x Z2 
invariance of Ui is also necessary if M.vf is to exhibit the TBM structure. This follows in 
a straightforward manner. Assume that M.uf has TBM structure of Eq. ([T]). The Upmns 
matrix in this case can be chosen to have the form in Eq. ([2|. Since Ui 
the following form in the basis specified by Eq 



U uU p jsf g 1 it has 



OtbaiPQ*OJ.^]^,j , 



^ 2pi +P2 P2- Pi P2 - Pi 

P2-P1 \ {Pi + 2p2 + 3^3) \ (pi + 2p2 - 3p3 
\ P2-P1 ^(Pi + 2p2 - 3p3) ^{pi + 2p2 + Sps) J 



(13) 



where pi denote the elements of the diagonal phase matrix PQ*. Interestingly, the above 
Ui is obtained from the general TBM Aiuf Eq. ([T]), by replacing the neutrino masses with 
the phases Pi and like Ai^j such a Ui is automatically Z2 x Z2 symmetric. Thus Eq. ^ 
also becomes necessary for the Z2 x Z2 invariance of Aiuf ■ Eq. (13) provides an alternative 



parametrization of f/;. It reduces to the earlier parametrization in Eq. ( 12 ) with the definition 



P3 = 



(14) 



with COST] = —3cgcosf3 and sin 77 = — cgsin/?. 

We note that the and Ui in Eqs. ( |8pl| are defined up to a simultaneous redefinition 
Uu — UUu and Ui — )■ UUi. In addition, Ui can also be multiplied by an arbitrary phase 
matrix from the right. Since Pi is arbitrary, the Ui in some model may not appear to have 
the Z2 X Z2 invariance even in a basis with chosen as in Eq. ([s]). But the above exercise 
shows that Ui can always be chosen to have the TBM form by appropriate rephasing of the 
charged lepton fields. 
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The above reasoning can be applied to more general patterns of mixing and not just to 
TBM. The key role in this construction is played by the fact that one can always choose 
a basis in which Mj, is Z2 x Z2 symmetric. This follows from the fact that the Z2 x Z2 
symmetry does not put any restrictions on the neutrino masses but only on the structure of 
mixing. As long as the neutrino mass matrices obey such "mass-independent" symmetries, 
the above construction of determining the most general Vi can be carried through. One can 
indeed define [TU] an appropriate Z2 x Z2 symmetry corresponding to every mixing pattern 
and then impose this symmetry on t/; to obtain the desired mixing structure in the flavour 
basis. As an example, in case of the /x-r symmetry, one can always choose a basis in which 
My is /i-r symmetric. Then requiring that Vi also be yU-r symmetric, we arrive at a general 
forms for Ui and which lead to a ^-t symmetric hAv] '■ 



f X A a\ 

ABC 
A C B 



U, 



(15) 



where Pil,r = diag.(l, e*^-^-^, e*^-^--^) and Ui has the same form of Eq. (12) but now 9 is an 
independent parameter and is not a function of (3 as before. 



III. 50(10) MODEL AND TBM 

We now integrate the above leptonic structures into an 5*0(10) model. We make the 
following assumptions which lead to simplification and allows us to obtain quantitative 
description. We (1) consider a supersymmetric S'O(IO) model with the Higgs transforming 
as 10, 126, 120 representations of 5*0(10) (2) impose the generalized parity as in pTI - ITS] 
leading to Hermitian mass matrices and (3) assume that the dominant contribution to M^, 
is a type-II seesaw, i.e. linear in the 126 Yukawa coupling. The last assumption with its 
attractive consequences is made in a number of 5*0(10) models [HI [121 HSl [T7]. The type-II 
dominance may be achieved by adding more Higgs fields. The type-I contribution can be 
suppressed by pushing the B — L breaking scale high while the type -II contribution can be 
relatively enhanced if the Higgs triplet remains below the GUT scale. The latter spoils the 
gauge coupling unification. It was pointed out in [18] that breaking 5*0(10) to SU{b) x f/(l) 
around 10^^ GeV can suppress the type-I contribution and presence of 54 plet can allow a 
complete multiplet transforming as 15 in 5*f/(5) to remain light. One can achieve in this 
way the type-II dominance without sacrificing the gauge coupling unification. The same 
phenomena was analyzed subsequently [IH] in a model with fields transforming as 54 + 45 
(instead of the standard 210) and having a 120 plet of Higgs. The fields 54 -|- 45 required to 
achieve the type-II dominance do not effect the fermion masses and the fermion sector can 
be the same as considered here. 

The fermion mass relations in this case after electroweak symmetry breaking can be 
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written in their most general forms as p^HTS] : 



Ma = H + F + iG , 
M„ = r{H + sF + itu G ), 
Mi = H-3F+ itiG , 
Md = r{H-3sF + itD G ), 
Ml = tlF , 

Mr = r-^^F (16) 

where (G) H, F are real (anti) symmetric matrices, r, s,ti,tu,tD,rL,rR are dimensionless 
real parameters. The effective neutrino mass matrix for three light neutrinos resulting after 
the seesaw mechanism can be written as 

M, = TlF - trMdF-^MI = M'J + . (17) 

The first term proportional to F denotes type-II seesaw contribution. In the numerical 
analysis that follows, we shall assume that My is entirely given by this term and subsequently 
analyze the effect of a small type-I corrections on the numerical solution found. 

We can always rotate the 16-plet fermions in generation space in such a way that oc F 
is diagonalized by the TBM matrix. 

F ^ R^FR = Ftbm = Otbm Diag.(/i, /a, h) O^^m (18) 

where fi are now real eigenvalues of F and the Otbm is given by Eq. (|3]). The matrix 
{G)H maintains its (anti)symmetric form in such basis and we use the same label for them 
in the rotated basis. In case of type-II seesaw dominance, the light neutrino mass matrix 
Mj^ = tlFtbm has the form given on the RHS of Eq. ([T]). The model has altogether 
17 independent real parameters (3 in Ftbm, 6 in iJ, 3 in G, r, s, t^, U and r^) which 
determine the entire 22 low energy observables of the fermion mass spectrum. Some of these 
parameters can be fixed by the known values of observables directly. As noted in [15], the 



5*0(10) relation for the charged lepton mass matrix in Eq. (16) can be rewritten as 

H + ttiG = ViDiV^ + 3FTBM , (19) 

where Di is a diagonal charged lepton mass matrix. Since H and G are real, the real and 
imaginary parts of the RHS separately determine H and tiG in terms of the charged lepton 
masses, parameters of Ftbm and Vi. V; is a unitary matrix that diagonalizes Mi and contains 
nine free parameters in the most general situation. One can suitably write Vi = ViP where 



P is diagonal phase matrix and Vi contains six real parameters. From Eq. (19), it is easy 
to see that the phase matrix P does not play any role in determining H and G and can be 
removed. So the nine real parameters of LHS can be related to six real parameters of Vi, 



three charged lepton masses and parameters of Ftbm in Eq. (19). This fixing helps us in 
numerical analysis as we will see in the next subsection. 
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We shall present numerical analysis in two different cases. (A) Corresponding to the most 
general Vi (B) with Vi = Ui given as in Eqs. ( [ll|l2| ). The case (A) has already been studied 
numerically in [111 1121 US] ■ We refine this analysis using a different numerical procedure and 
taking into account the results of the most recent global fits [S] to neutrino data. This also 
serves as a benchmark with which to compare the case (B) which leads to the exact TBM 
at Mgut- 



A. Numerical Analysis: The most general case 



We study the viability of Eq. (16) with the experimentally observed values of fermion 
masses and mixing angles through numerical analysis. For this, we construct a function 
defined as 

'P^-Oi ^ 



X 



cr,; 



(20) 



where the sum runs over different observables. P,- denote the theoretical values of observ- 



ables determined by the expressions given in Eq. (16) and Oi are the experimental values 
extrapolated to the GUT scale, denote the la errors in Oj. Our choice of the input 
values of quark and lepton masses and quark mixing angles are the same as used in [11] . In 
this data set, the charged fermion masses at the GUT scale are [20] obtained from the low 
energy values using MSSM and tan/3 = 10. We use the input values of lepton mixing angles 
from [3] which includes results from T2K and MINOS. The effect of the RG evolution on 
the quark mixing angles is known to be negligible. This is also true for the lepton mixing 
angles in case of the hierarchical neutrino spectrum. We assume such hierarchy in neutrino 
masses and therefore the input values of the quark mixing angles, CP phase and neutrino 
parameters we use correspond to their values at low energy. We reproduce all these input 
values in Table [T] for convenience of the reader. 



GUT scale values with propagated uncertainty 


md(MeV) 


1.24 ±0.41 


AmL(eV2) 


(7.58 ±0.22) X 10"^ 


ms(MeV) 


21.7 ±5.2 


Am2,^(eV2) 


(2.35 ±0.12) X 10~3 


m;,(GeV) 




sin 6*^2 


0.2243 ±0.0016 


m„(MeV) 


0.55 ±0.25 


sin 6*^3 


0.0351 ±0.0013 


mc(GeV) 


0.210 ±0.021 


sin 0^3 


0.0032 ± 0.0005 


mt(GeV) 


82.4130-3 


sin2 


LI.OUD_o 015 


me(MeV) 


0.3585 ± 0.0003 


s\v? 6*^3 




m^(MeV) 


75.672 ± 0.058 


sin^ ^^3 


021+°-°°'^ 

U.U^i_Q QQg 


m^(GeV) 


1.2922 ±0.0013 


Jcp 


(2.2 ±0.6) X 10-5 



TABLE I. Input values for quark and leptonic masses and mixing angles in the MSSM extrapolated 
at Mgut = 2 x 10^^ GeV for tan /3 = 10 which we use in our numerical analysis. 
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We fit the above data to the fermion mass relations (16) predicted in the model by 



numerically minimizing the function. As already mentioned above, this exercise has been 
done in jlT] recently and a very good fit corresponding to Xmm = 0.127 is found. We repeat 
the same analysis because of the following differences in our fitting procedure and because 
of the emergence of new results on 6*^3 [5], E]. 



Compared to other observables, the charged lepton masses are known very precisely 
with extremely small errors in their measurements. Instead of fitting them through 



minimization, we use their central values as inputs on the RHS of Eq. ( 19 ). Because of 



this, our definition of the function in Eq. (20) does not include the charged lepton 
masses in it. 



• We also use the central value of the solar to the atmospheric mass squared difference 
ratio as an input and use it to fix /a through the following relation: 

(2\ 

After obtaining the solution, the overall scale of neutrino masses r^, at the minimum 
is determined by using the atmospheric scale as a normalization. 

As a result of these simplifications, the function in our approach includes only 13 observ- 
ables, 6 quark masses, 3 quark mixing angles, a CKM phase and 3 lepton mixing angles. 
These are complex nonlinear functions of 12 real parameters {2 in Ftbm, 6 in V;, r, s, ti and 
tu)- This is numerically minimized using the function minimization tool MINUIT. The 
results of our analysis are shown in column A in Table [Tij 

We obtain an excellent fit corresponding to Xmin = 0.0061 which is significantly better 
than Xmin = 0.127 obtained in [11] using different procedure and different data set. Pa- 
rameters obtained for the best fit solution are shown in Appendix. All the observables are 
fitted within the O.OScr deviation from their central values. The solution at its minimum 
almost reproduces the central value of sin^ 6*^3 obtained in the global fit [3]. We also de- 
termine a parameter dpT introduced in [H] which quantitatively measures the amount of 
fine-tuning needed in the parameters for obtained fit when compared with doata which is a 
similar parameter obtained from the data only. From our fit, we obtain dpT ~ 2.9 x 10^ 
compared to dijata ~ 220. These parameters depend on the definition of and since we 
do not include the charged lepton masses -which have very small errors- in our x^, both the 
parameters dpr and doata are an order of magnitude smaller than jTl]. However the ratio 
dData/dpT ~ 7.6 X 10~^ obtained from our fit is almost similar to dnataldpT ~ 8 x 10~^ 
obtained in [TT] which shows that both of the solutions need substantial level of fine tunings 
in the model parameters. 



8 





Case A 


Case Bl 


Case B2 


Observables 


Fitted value 


Pull 


Fitted value 


Pull 


Fitted value 


Pull 


mdiMeVl 


1.2339 


-0.0148738 


1.22098 


-0.0463899 


1.02686 


-0.519852 


mJMeVl 

L J 


21.7214 


0.00411949 


21.9922 


0.0561874 


22.0058 


0.058806 


mJGeVl 

u L J 


1.06614 


0.0438763 


1.16345 


0.738942 


1.2842 


1.60145 


m„[MeV] 




000073755 


550234 


n 000936368 


n 550787 


0031 4771 


mJGeV] 


0.209977 


-0.00111886 


0.209952 


-0.00230315 


0.210481 


0.0229054 


mt[GeY] 


OZ.OZYO 


(J.UU42i748 


82.5855 


U. 1)1)512198 


81.7487 


— U.(J44(JUoz 


me [Me V J 


0.3585 


— 


0.3585 


— 


0.3585 


— 


m^[MeV] 


75.672 


— 


75.672 


— 


75.672 


— 


m^[GeV] 


1.2922 


— 


1.2922 


— 


1.2922 


— 




0.0323 


- 


0.031875 


- 


0.031875 


- 


sin 


0.224299 


-0.000688878 


0.2243 


0.0002182 


0.224303 


0.0019076 


sin ^23 


0.0351032 


0.00246952 


0.0350951 


-0.0038047 


0.0351294 


0.022597 


sin 


0.00320513 


0.0102511 


0.00319436 


-0.0112796 


0.0031749 


-0.0502087 


sin^ 9 


0.306119 


0.00660722 


0.3333 




0.3333 




sin^ OL 

Zo 


0.418475 


-0.0508353 


0.5 




0.5 




sin^ 05_3 


0.0207708 


-0.0286467 












Jcp 


2.19 X 10-^ 


-0.0183401 


2.21 X 10"^ 


0.0194165 


2.25 X 10"^ 


0.0845729 


^MNS 


282.396 

















ai 


142.975 


— 


160.829 


— 


180 


— 


Q!2 


22.4851 




318.593 











6.89 X 10~i° 




9.82 X 10-1° 




3.53 X 10-9 




■ 




0.0061 




0.5519 




2.8510 


dpT 




28751 




88546 




31171 


doata 




220 




200 




200 



TABLE II. Best fit solutions for fermion masses and mixing obtained assuming type-II seesaw 
dominance in the SUSY 50(10) model with 10 + 126 + 120 Higgs. Various observables and their 
pulls at the minimum are shown for three different cases correspond to (A) the general (non TBM) 
leptonic mixing, (Bl) Exact TBM leptonic mixing with Ui of Eqs. ( |11|12| ) and (B2) Exact TBM 
leptonic mixing with diagonal Mi and ti = (See the discussions in the text for more details). The 
predictions of different approaches are shown in boldface. 



B. Numerical Analysis: Exact TBM 



After discussion of the above general case, we now specialize to the case of the exact 
TBM. This case is of considerable theoretical interest since it can point to some underlying 
symmetry existing at Mqut- We can implement the exact TBM in a model independent way 
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by choosing Vi = Ui in Eq. (19). With this choice, all the leptonic mixing angles get fixed to 
their TBM values. Also the central value of the ratio of the solar to the atmospheric (mass)^ 
differences is used as input and a parameter tl is determined at the minimum by using the 
atmospheric scale. Thus the function in Eq. (20) now involves only observables in the 
quark sector. As already discussed, Ui in Eq. ( ll|l2 ) is parameterized by three phase angles 
a, (3 and 5. An overall phase a is irrelevant for the physical observables and can be removed. 
This leaves only 8 real parameters (2 in Ftbm, 2 in Ui, r, s, ti and ty) which are fitted 
to the 10 observables in the quark sector by minimizing the x^- The results are shown in 
column Bl in Table [lH The obtained fit corresponds to Xmin — 0.552 (x^^„/d.o.f. = 0.276). 
Only the fitted value of mb deviates slightly from the central value with a 0.74a pull. All 
the remaining observables are fitted within 0.06a. A set of parameters obtained for this 
solutions are shown in Appendix. The fit obtained here is not significantly different from 
the general case discussed before showing that all the fermion masses and mixing angles can 
be nicely reproduced along with the exact TBM within the S'0( 10) framework discussed 
here. 

Before we discuss possible perturbations in TBM pattern, let us discuss a very special 
case corresponding to a diagonal Mi. This corresponds to Ui coinciding with an identity 
matrix and is a special case of Eq. ^Vj^ with /3 = f and 5 = 0. Since Mi is real and 
diagonal in this case, tiG must vanish in Eq. (19). If G = then the quark mass matrices 
also become real and there is no room for CP violation. The viable scenario must therefore 
have nonzero G and hence ti = 0. As a result, unlike before, the three parameters in G do 
not get determined from Mi, see, Eq. (19) and remain free. They can be fitted from the 
quark sector observables. We carried out a separate numerical analysis for this particular 
case and the results are shown in column B2 in Table |TT} The fit obtained gives relatively 
large Xmin = 2.85 (Xmin/d-o.f. = 1.43) with more than la deviation in the bottom quark 
mass. Although the obtained Xmin statistically acceptable at 90% confidence level, it is 
not as good as the previous one and we shall not consider this case with Ui = I any further. 



IV. PERTURBED TBM 

The TBM is an ideal situation and various perturbations to this can arise in the model. 
We need to analyze these perturbations in order to distinguish this case from the generic 
case without the built in TBM. A deviation from tri-bimaximality can arise due to 

1. renormalization group evolution (RGE) from Mqut to Mz. 



2. small contribution from the sub dominant type-I seesaw term in Eq. (17). 

3. the breaking of the Z2 x Z2 symmetry in Ui which ensured TBM. 

The effect of (1) is known to be negligible ^21] in case of the hierarchical neutrino mass 
spectrum which we obtain here. We quantitatively discuss the implications of the other two 
scenarios via detailed numerical analysis in the following subsections. 
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A. Perturbation from type-I seesaw 



Depending on the GUT symmetry breaking pattern and parameters in the superpotential 
of the theory, a type-I seesaw contribution can be dominant or sub dominant compared to 
type-II but it is always present and can generate deviations in an exact TBM mixing pattern 
in general. In the approach pursued here it is assumed that such contribution remains sub 



dominant and generates a small perturbation in dominant type-II spectrum. Eq. (17) can 
be rewritten as 

M, = rL{F-^MDF-^Ml) (22) 

where ^ = tr/tl determine the relative contribution of type-I term in the neutrino mass 
matrix. 



The second term in Eq. (17) brings in two new parameters ^ and to present in the 
definition of M/) in Eq. (16). These parameters however affect only the neutrino sector. We 
isolate the effect of type-I contribution by choosing other parameters at the minimum 
found in Section III-B. ^ and to remain unconstrained at this minimum and their values 
do not change the obtained earlier since the latter contains only the observables in the 
quark sector. The ^, tz? however generate departure from the exact TBM. We randomly vary 
the parameters and t/? and evaluate the neutrino masses and mixing angles. While doing 
this, we take care that all these observables remain within their present Sa [3] limits. Such 
constrains allow very small values of |^| < 10~^. The correlations between different leptonic 
mixing angles found from such analysis are shown in Fig. ([T|. 

It is seen from Fig. ([T]) that the perturbation induced by type-I term cannot generate 
considerable deviation in the reactor angle if the other two mixing angles are to remain 
within their Sex range. In particular, requiring that sin^ 6*^2 remains within the So" range 
puts an upper bound sin^^^g < 0.0002 which does not agree with the latest results from 
T2K and MINOS showing that a small perturbation from type-I term cannot be consistent 
with data when the type-II term displays exact TBM. 



B. Perturbation from the charged lepton mixing 

A different class of perturbation to TBM arise when Ui deviates from its Z2 x Z2 symmetric 



form given in Eq. (11). In this case, the neutrino mass matrix has TBM structure but 
the charged lepton mixing leads to departure from it. This case has been considered in 
the general context [22] as well as in 5*0(10) context [TT]. Within our approach, we can 
systematically look at the perturbations which change the values of any one or more angles 



from the TBM value. For example, we can choose Ui as given in Eq. (15) and look at the 
quality of fits in this case compared to the exact TBM solution. This choice (corresponding 
to ii-T symmetry) leaves ^23 cind 6*^3 unchanged but perturbs ^[2- Alternative possibility is 
to simultaneously perturb all three mixing angles and look at the quality of fit compared to 
the exact TBM case. We follow this approach. For this we choose Ui to be a general unitary 
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0.35 0.40 0.45 0.50 0.55 0.60 

sin^eL 




FIG. 1. Correlations among the lepton mixing angles when two real parameters ^ and t£) are 
varied randomly. The points with different colors correspond to {tol < 1 (blue), 1 < I^dI < 5 (red) 
and 5 < < 10 (green). The black lines are the updated central values of sin^ ^-i^d sin^ 
obtained by the global fits on neutrino oscillation data [3j. The unshaded and the shaded regions 
correspond to la and 3a bounds respectively. 



matrix and go back to the analysis in Section III-A. There, we have fitted the solar and the 
atmospheric mixing angles to their low energy values given in Table |T) Here, we modify the 
definition of pin down a specific value Oi of the mixing angles pi by adding a term 



to Xg that contains all the observables of the quark sector. Sum in Eq. (23) runs over the 
three lepton mixing angles. The = Xg + Xim then numerically minimized to fit the 
13 observables determined in terms of 12 real parameters as mentioned in Section III-A. 



Artificially introduced small errors in Eq. (23) fix the value Oj for pi at the minimum of the 
X^. We then look at the quantity 

(24) 

which represents the fit to the quark spectrum when the lepton mixing angles Pi are pinned 
down to values Oj. We repeat such analysis by randomly varying Oi within the allowed 3cr 
ranges of lepton mixing angles j3]. The results are displayed in Fig. ([2]). 
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FIG. 2. Correlations among the lepton mixing angles in case of the most general charged lepton 
mixing matrix Ui. The points with different colors correspond to Xmin < (green), 1 < Xmin 
(blue) and Xmin — ^ (I'ed). The black lines are the updated central values of sin^ ^'^^g S'^d sin^ 
obtained by the global fits on neutrino oscillation data |3j. The unshaded and the shaded regions 
correspond to la and 3a bounds respectively. 



We plot the correlations among the lepton mixing angles and show the corresponding 
values of Xmin three different regions. The points corresponding to Xmin < 1 (green) 
represent very good fit in which all the observables are fitted within la. The obtained fit 
shown by the points corresponding to 1 < Xmin < ^ (blue) is not as good as the previous one 
but it is statistically acceptable. The points for xLin > 4 (red) represent poor fit and can 
be ruled out at 95% confidence level. Fig. ^ shows definite correlations between 6*23 and 
9[2. It is seen that the region Xmin < 4 falls largely below sin^ O^-^ = 0.58 for sin^ = 0.306. 
It is also seen from the figure that the entire range 0.001 < sin^ < 0.044 is consistent 
with statistically acceptable fits to fermion spectrum. This is to be contrasted with the 
previous case where perturbation from type-I seesaw term led to an upper bound. The 
bounds obtained numerically allows us to clearly distinguish the case of the exact TBM at 
Mgut in comparison to the one in which the charged leptons lead to departures from the 
tri-bimaximality. 
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V. SUMMARY 



The presently available information on the leptonic mixing may be described by a TBM 
structure for Aiuf combined with appropriate perturbations generating relatively large Ol^. 
We have analyzed the viability of this scenario in a larger context of the grand unified 
SO{10) theory taking a specific model as an example. The TBM structure for the neutrino 
mass matrix is a matter of choice of the basis [H]. Thus the existence of TBM is linked to 
the structure of the charged lepton mixing matrix Ui in this basis. All the related studies in 
this context [HI [22] assume that Ui deviates slightly from identity and discuss the breaking 
of the TBM pattern through such Ui. We have shown that it is possible to construct a class 
of non-trivial Ui quite different from identity which preserve the TBM structure of M^, when 
transformed to the flavour basis. Identification of such non-trivial Ui becomes crucial in the 
context of 5*0(10) and allows us to obtain a viable fit to fermion spectrum keeping TBM 
intact. The quality of fit obtained in this case is excellent as shown in Table |TT] and differs 
only marginally from a general situation without imposing the TBM structure at the outset. 

The existence of TBM at the GUT scale may be inferred by considering its breaking which 
can arise in the model and the reactor mixing angle is a good pointer to this. The quantum 
corrections are known ^I] to lead to very small for the hierarchical neutrinos. Simi- 
larly, corrections coming from type-I seesaw term imply an upper bound, sin^6'^3 < 0.0002 
as discussed in Section IV-B. These two cases can be ruled out by relatively large value 
of as indicated by the observations from T2K and MINOS. These cases are in sharp 
contrast to a situation in which one does not impose the TBM at Mqut and determined 
Ui from a detailed fits to fermion masses. We have analyzed this scenario using the results 
from the latest [3] global fits to neutrino oscillation data. It is found that the entire range 
0.001 < sin^ < 0.044 is consistent with the detailed description of all the fermion masses 
and mixing angles. 
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VI. APPENDIX 

In this Appendix, we show the parameter values for each of the two cases (A) and (Bl) 
in Table [Tij The given values of the matrices H, F, G and parameters r, s, ti and t„ are at 
the minimum of and obtained from our fitting procedure. All the physical mass matrices 



can be constructed from the given numbers using the 5*0(10) relations in Eq. (16). 
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Case A: The most general Vi 



The parameter values corresponding to the best fit solution shown in column A in Table 
[II] are the following. 

/ 0.00286198 0.00178975 -0.0399313 
H = 0.00178975 0.0196685 0.00741417 I GeV 
-0.0399313 0.00741417 1.10192 



G 



-0.00235758 
-0.0053869 
-0.0053869 





-0.0053869 
-0.0394516 
0.0317071 



-0.00144676 
0.00144676 
V 0.00442083 0.0183052 



-0.0053869 
0.0317071 
-0.0394516 

-0.00442083 
-0.0183052 




GeV 



GeV 



(r, s, tu, ti) = (76.2076, 0.54487, - 0.870094, - 14.5065) 



Case B: Exact TBM lepton mixing (Vj = Ui) 



The parameter values corresponding to the best fit solution shown in column Bl in Table 
[IT] are the following. 

/ 0.00111546 0.000428794 -0.020097 
H = 0.000428794 0.0188148 0.0639684 I GeV 
-0.020097 0.0639684 1.17731 



V 
/ 

V 



-0.0016484 
-0.00375321 
-0.00375321 



G 



( ° 

0.00278749 

\ 0.0378996 



-0.00375321 
-0.0276748 
0.0222732 

-0.00278749 - 


0.0820076 



-0.00375321 
0.0222732 
-0.0276748 

-0.0378996 \ 
-0.0820076 




GeV 



GeV 



(r, s, tu, ti) = (70.5742, 0.526782, 0.551405, 2.15986) 
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